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ABSTRACT 

The purpose of this research was to compare numerical predictions with experimental 
results and to devise an accurate laser anemometry technique to measure in the endwall 
region of a confined annulus. Flowfield characteristics were computed using a three- 
dimensional flow solver with the numerical plane coincident with the experimental 
measurement location. An annular turbine cascade, designed for laser-Doppler velocimetry, 
was modified to obtain blade passage midspan surface pressure measurements. A range of 
predicted subsonic and transonic midspan surface pressures were compared favorably with 
experimental measurements. Two-dimensional Mach number, flow angle, and turbulence 
intensity measurements were obtained with a fiber-optics laser-Doppler velocimeter. The 
measurements were performed through a 1.0922 millimeter opening in the endwall at depths 
ranging from 0.01 mm to 3.34 mm and the results were compared with numerical 
predictions. 
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I. INTRODUCTION 


Advanced aeropropulsion systems have required substantial technological 
improvements in turbomachinery. These depend on the availability of complex analysis and 
experimental tools. A thorough understanding of the generation and development of 
secondary flows in annular blade rows has provided insight into the design and performance 
of turbomachinery. Secondary flows, together with tip leakage flows, produce considerable 
flow distortions and losses in the endwall region. [Ref. 1] These losses can be minimized and 
turbomachinery efficiency can be improved with a more precise understanding of the flow 
mechanics and the ability to numerically predict the flow field. The current emphasis on 
turbomachinery design centers on numerical analysis. 

This report specifically documents the investigation of the flow through an annular 
turbine cascade (ATC). References 2 through 4 include continuing research in the field of 
laser-Doppler velocimetry (LDV) and numerical prediction in a large annular turbine 
cascade. The authors provided experimental two and three-dimensional velocity, flow angle, 
and turbulence intensity at constant axial positions within an ATC passage. These 
measurements, in addition to vane surface static pressure measurements, supplied a test case 
for three-dimensional turbomachinery computer programs. They concentrated on obtaining 
comparisons ahead of, inside, and downstream of the blade passage and concluded that the 
largest difference between experimental and computational results was in the endwall region 
where viscous and secondary flow effects were the greatest. 

Two previous investigations of the flow through the ATC have been conducted at the 
Naval Postgraduate School. Reference 5 included design and manufacturing information of 
the annular turbine cascade that was developed to determine the limitations of LDV 
measurements in a confined annulus. Reference 6 included additional laser and pressure 
probe access modifications and initial LDV measurements to the same ATC. This report 
includes further ATC modifications for midspan blade surface pressure measurements and 
LDV measurement techniques of the endwall flow. Radial two-dimensional fiber-optic probe 
traverses were performed, through a small access hole in the outer casing, to coincide with the 
numerical exit plane. Circumferential surveys were obtained at different radial locations close 
to the endwall. Blade midspan surface pressures were measured within one blade passage at 
various inlet total-to-downstream hub-static pressure ratios. Blade surface pressure and 
endwall flow measurements were compared with numerical predictions obtained using a 
three-dimensional viscous computer program. 



The two most noteworthy comparisons were at a subsonic pressure ratio and a 
pressure ratio corresponding to sonic exit conditions. Comparisons with the LDV data were 
performed at the subsonic flow condition. The numerical blade surface pressure distributions 
compared well with the experimental results, particularly for the sonic exit condition for 
which trailing edge shocks were predicted. 
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H. EXPERIMENTAL SETUP 



A. TEST FACILITY AND ANNULAR TURBINE CASCADE 

Airflow for the annular turbine cascade experiment was provided by a VA-312 Allis- 
Chalmers 12-stage axial-flow compressor located at the Turbopropulsion Laboratory of the 
Naval Postgraduate School. The compressor was operated at 12,000 rpm at various discharge 
pressures and provided metered air to a plenum chamber. Air from the plenum was routed to 
a 232.918 mm (9.170 in) diameter bellmouth and test section through honeycomb flow 
straighteners in a 254 mm (10 in) flanged steel pipe as shown in Figures 1 and 2. 


Figure 1. Front View of Experimental Apparatus 

Flow stagnation pressure was measured at two upstream locations. One combination 
probe provided pressure-setting information to a mercury manometer board and a digital 
readout of flow stagnation temperature, while the second probe was connected to a Scanivalve 
(Figure 2, lower right). Four (averaged) upstream static ports and four (averaged) inner hub 









downstream (one-half axial chord) static ports were also connected to a Scanivalve for 
automated pressure data acquisition. 


Figure 2. Top View of Experimental Apparatus 

Atomized glycerin particles constituted the LDV seed material which were introduced 
through a 7.938 mm (0.313 in) diameter copper tube approximately ICS tube diameters 
upstream of the test section. Seed atomization was performed using a commercial TS1, Inc., 
Six-Jet Atomizer connected to the laboratory compressed air supply. Laser traverse 
mechanism buffet, caused by the exiting flow, was minimized with a prefabricated aluminum 
cone (Figure 2, upper center) and Reference 6 showed that flow characteristics were unaltered 
with the cone attached. 

Each blade was designed with a combination of simple circular arcs and line segments 
and included a leading edge radius of 2.413 mm (0.095 in), trailing edge radius of 0.305 
mm (0.012 in), and axial chord of 24.77 mm (0.975 in). The annular stator row was 
manufactured from 2218-T61 aluminum and consisted of 31 blades with a midspan spacing 



of 21.77 mm (0.857 in), resulting in a blade solidity of 1.14. The inner hub radius was 
98.93 mm (3.895 in) and the outer case radius was 116.46 mm (4.585 in) with the same 
profile at all radii. Reference 5 included the original set of manufacturing drawings and 
Reference 6 included a description of the wake positioning system. 

B. PRESSURE MEASUREMENT MODIFICATIONS 

The original ATC did not provide the capability to measure blade surface pressures. 
Within a single passage, seven suction-side static ports and four pressure-side static ports were 
drilled orthogonal to the blade surface at midspan. Figures 3 and 4 show front and rear views 
of the instrumented blades. As shown in Figure 5, each static port was 0.406 mm (0.016 in) 
in diameter and each spanwise hole was 1.321 mm (0.052 in) in diameter. 

All spanwise holes were sealed at the tips and stainless steel tubes were cemented into 
the hub openings. The tubes were connected to the Scanivalve with plastic tubing which was 
fed out through a sting from the center body. (Figures 1 and 2) Each port's circumferential 
position was measured using a 2.375 mm (0.0935 in) diameter stylus and later converted 
graphically to an axial chord position to allow computational comparisons. Figure 5 shows 
port numbering and dimensions. 



Figure 3. Blade Leading Edge View With Pressure Measurement Modifications 




Figure 4. Blade Trailing Edge View With Pressure Measurement Modifications 



Figure 5. Blade Static Port Numbering Sequence And Dimensions 







C. PRESSURE DATA ACQUISITION 

The data acquisition system, for the pressure measurements, is shown schematically in 
Figure 6. All data acquisition was remotely controlled by a Hewlett-Packard 9000 computer 
system. Appendix A contains the program utilized to conduct all pressure data acquisition. 

A Scanivalve was connected to a Model HG-78K Scanivalve controller, which in turn was 
connected to a Hewlett-Packard, Model 3456A Digital Voltmeter and Model 3495A Scanner 
via a HP-IB instrument bus. Scanivalve calibration was performed to within an accuracy of 
+/- 0.1 inches mercury. Table A1 in Appendix A relates each Scanivalve port to its respective 
pressure measurement. 



Figure 6. Pressure Data Acquisition Schematic 










































D. LASER-DOPPLER VELOCIMETER 

Part of the laser apparatus is shown in Figure 7. The probe, processor, computer, and 
traverse mechanism are shown in Figures 1 and 2. The LDV system is shown schematically in 
Figure 8. A LEXEL Model 95 four-Watt argon-ion laser was connected to a TSI, Inc., Model 
9201 ColorBurst multicolor beam separator. The beam separator divided the incoming light 
into shifted and unshifted beams, with the shifted beam receiving a 40 MHz frequency shift 
from a Bragg cell. The two beams were further split into three polarized pairs: green (514.5 
nm), blue (488 nm), and violet (476.5 nm). 

Individual couplers on the ColorBurst directed each beam to the laser probe via a 
fiber-optic cable. Each fiber-optic probe contained receiving optics which directed the return 
signal to a TSI, Inc., Model 9230 ColorLink multicolor receiver. The ColorLink provided 
photomultiplier and frequency-shifting functions. All conditioned ColorLink signals were 
sent to a TSI, Inc., IFA-750 digital burst correlator where valid Doppler signals were 
identified and digitized. 

The fiber-optic probes were mounted to a LINTECH, Model 41583 traverse table. 

An Applied Motion Products System 1618 traverse controller was used manually to control 
traverse table movement. All ColorLink and IFA-750 functions and LDV data processing 
were accomplished remotely by computer using TSI's menu-driven software, FIND (FLOW 
INFORMATION DISPLAY) version 4.04. 

E. EXPERIMENTAL PROCEDURE 

1. Midspan Surface Pressure Measurements 

Midspan surface pressure measurements were obtained with the pressure data 
acquisition system. The pressure ratio (Prat) was defined as the downstream hub-static 
pressure (Phub) divided by the upstream stagnation pressure (Po)- Each pressure ratio was set 
by metering the upstream stagnation pressure until a desired mercury manometer column 
height was achieved. Five pressure ratios (0.5070, 0.6041, 0.6815, 0.8077, and 0.9054) were 
considered and during each run all the blade surface pressures, Po, and Phub were recorded. 

2. Laser Alignment 

LDV alignment for endwall flow measurements was accomplished as shown in Figure 
9. The objective of the LDV alignment procedure was to center the probe volume in the 
1.0922 mm optical access hole at a known and repeatable radial distance. The two- 
dimensional fiber-optic probe was attached to a mounting-bracket micrometer which allowed 




Figure 7. LDV Bread Board (with Laser, Color Separator, and Receiving Optics Module) 



Figure 8. LDV System Schematic 



















































probe radial travel in 0.01 mm increments. After zeroing the micrometer, the traverse table 
was manually advanced forward until the four beam separation was minimized, yet discernible 
with the naked eye (Figure 9, sketch A). The digital position-indicator reading on the 
traverse controller was noted and the process repeated for forward travel until the same image 
reappeared. The midpoint of the traverse table positions was defined as the center of the 
probe volume with the face of the optical access plate as a radial reference point. Reference 6 
described complete dimensions of the optical access plate. 



Optical Access Hole 

1 1.0922 mm diameter _ 


1 


: t 

3 

1 A 

B 


Figure 9. LDV Alignment Schematic 

Horizontal and vertical alignment positions (Figure 9, sketch B) were obtained by 
noting the digital position indicator on the traverse controller as the probe volume touched 
the left(l), right(2), bottom(3), and top(4) inner edge of the optical access hole. The 
horizontal and vertical reference positions were defined as the center of the optical access 
hole. The probe was then traversed inward by 1.14 mm (the thickness of the optical access 
plate), at which point the center of the probe volume coincided with the outer (case) wall of 
the turbine cascade. This alignment technique was repeatable and ensured that the probe 
volume passed cleanly through the center of the optical access hole. 

3. Endwall Measurements 

All endwall measurements were conducted at a pressure ratio of 0.9054. The laser 
beams were aligned with the downstream optical access hole as described above, and then 
traversed manually using the mounting-bracket micrometer. Radial endwall surveys at 0.01, 
0.06, 0.18, 0.42, 0.89, 1.78, and 3.34 millimeters from the tip casing were conducted for 
peripheral (wake) angular settings ranging from -8 degrees to +8 degrees. One degree 
(wake) position increments were achieved by circumferentially rotating the blade row and 
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center body within the outer casing. At each circumferential position the inner section was 
secured in place with a locking bolt arrangement, the design of which was documented in 
Reference 6. Seeding position, ColorLink, and IFA-750 settings were adjusted for an 
optimum LDV data rate and minimum noise. All settings ensured a minimum data rate of 
100 samples per second, however; LDV data rate was extremely sensitive to position of the 
wand which introduced seed into the flow. 

A 'random' mode processor setting allowed a total of 1,024 samples between the green 
and blue channels with no user control over sample distribution. The 'coincidence' mode 
acquired 1,024 samples for each channel. Repeatability measurements for endwall velocity, 
flow angle, and turbulence intensity were separately conducted at three wake positions for 
both modes. A minimum pressure ratio was determined before data rate conditions became 
unacceptable. Raw data were converted within FIND and manually transferred to a 
spreadsheet for further processing. 
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m. COMPUTATIONAL FLUID DYNAMICS 


A. GRID GENERATION 

Grid generation was completed using the FORTRAN language program 
"Turbomachinery C GRID (TCGRID)" [Ref. 7], TCGRID inputs consisted of four lines of 
namelist inputs followed by a title, hub and tip geometry, and blade geometry. The blade 
geometry was input in cylindrical coordinates (z, theta, r) starting at the blade trailing edge 
and wrapping clockwise. The blade inputs were completed in slacked sections from the hub 
to the tip. The grid used, with resolution of 150 x 31 x 65 (i, j, k), was the same as that 
generated in Reference 6. The i-index was defined clockwise from the lower (pressure 
surface) exit to the upper (suction surface) exit, the j-index was defined from the blade 
surface to the periodic boundary, and the k-index was defined from the hub to the tip. All 
computational solutions were based on this gird. The final grid is shown in Figure 10. 
Appendix B contains the grid namelist input file and Figure B1 shows the blade geometry. 

B. COMPUTATIONAL SCHEME 

Flowfield solutions were obtained using two versions (206 and 208) of "Rotor Viscous 
Code 3-D (RVC3D)", a FORTRAN language program designed for analysis of three- 
dimensional viscous flows in turbomachinery. RVC3D was written to solve the thin-layer 
Navier-Stokes equations in Cartesian coordinates. The equations were discretized using 
second-order finite-differences in space and solved in time with a fourth-order Runge-Kutta 
scheme. Streamwise viscous terms were neglected using a thin-layer assumption, but cross¬ 
channel viscous terms were retained. A spatially-varying time step and implicit residual 
smoothing were used to accelerate convergence. [Refs. 8-10] Turbulence effects were 
modeled using a 3-D adaptation of the Baldwin-Lomax model and the Cebeci-Smith model. 

RVC3D version 206 was used to predict the flowfield for pressure ratios of 0.6041, 
0.6815, 0.8077, and 0.9054. The turbulence in the flowfield was computed with an 
adaptation of the Cebeci-Smith turbulence model. This version of RVC3D only 
accommodated subsonic exit boundary conditions whereby the hub static pressure was held at 
the pressure ratio and radial equilibrium was solved for the spanwise pressure distribution. 
Version 208 of RVC3D was used for the 0.5070 pressure ratio, which produced trailing edge 
shocks that extended to the exit plane. This newer version allowed for a supersonic exit 
boundary condition where the exit conditions were based on Giles' characteristic boundary 
conditions [Ref. 11]. The Baldwin-Lomax turbulence model was used for this test case since 


13 




Figure 10. Three-Dimensional Multiple Grid (150 x 31 x 65) 






a 'bug' was discovered that would not allow the code to properly recognize the Cebeci-Smith 
turbulence model [Ref. 12]. The code's author has since corrected the program. 

Appendix C contains an example namelist input file used to obtain a flow solution 
and a description of the steps required to run the code on the Naval Postgraduate School 
(NPS) and National Aerodynamic Simulation Program (NAS) Cray supercomputers. Output 
solution file (fort.3) information was visually examined with FAST and PLOT3D graphics 
software [Refs. 13 and 14]. Solution residual files (fort.4), blade surface pressure files 
(fort.7), and Mach number and flow angle files (fort. 7-11) were calculated from the solution 
file with the FORTRAN programs "pxy.F 1 and "plane.f" (Appendix G) for final graphical 
representation. 


15 



16 



IV. RESULTS AND DISCUSSION 


A. BLADE MIDSPAN SURFACE PRESSURE MEASUREMENTS 

Experimental blade-surface pressure measurements were averaged and a maximum 
repeatability error was calculated. Each blade static port reading was non-dimensionalized by 
the upstream stagnation pressure (Po). In order to provide an experimental comparison with 
Reference 6 for the 0.6815 pressure ratio, experimental repeatability was confirmed by 
conducting six runs; three runs each on separate days (Tables D1 and D2 of Appendix D). 
The maximum repeatability error was 0.87%. All other mean pressure data were based on a 
three run average at each pressure ratio and resulted in a maximum repeatability error 
ranging from 0.66% at the lowest pressure ratio down to 0.22% at the highest pressure ratio. 
Table 1 contains each midspan static port location and non-dimensional pressure 
measurement for all experimental pressure ratios. Tables D3 through D6 of Appendix D 
contain pressure data for pressure ratios of 0.5070, 0.6041, 0.8077, and 0.9054 respectively. 


Static Port 
.Number 

Fraction of 
Axial Chord 

Pressure 
Ratio = 

Pressure 
Ratio = 

Pressure 
Ratio = 

Pressure 
Ratio = 

Pressure 
Ratio = 
0.9(154 

1 

0.0000 


1.0019 

0.9999 

0.9996 


2 

0.1322 

0.9628 

0.9576 

0.9600 

0.9729 

0.9872 


0.4615 

0.8613 

0.8646 

0.8791 

0.9171 


4 

0.6604 

0.6986 

0.7223 

0.7595 

0.8433 

0.9195 

5 

0.7757 

0.5419 

0.6216 

0.6910 

0.8113 

0.9057 


0.8542 

0.5584 

0.6358 

0.7031 

0.8190 


7 

0.9183 

0.5200 

0.6166 

0.6928 

0.8169 




1.0082 

0.9987 





0.1322 

0.9955 

0.9878 

0.9893 

0.9923 

0.9973 

10 

0.6274 

0.9389 

0.9362 

0.9439 

0.9600 

0.9810 

11 

0.7508 

0.8853 

0.8874 

0.9015 

0.9327 

0.9668 


Table 1. Non-Dimensional Midspan Surface Pressure (P/Pq) 


B. LASER-DOPPLER VELOCIMETRY MEASUREMENTS 

Two-dimensional LDV measurements were obtained to an approximate depth of 3.34 
mm through a 1.0922 mm (0.043 in) diameter hole located one-half axial chord 
downstream. All LDV data were acquired at a pressure ratio of 0.9054 and are tabulated in 
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Tables El through E5 of Appendix E. All data resulted from programmed editing of 
histograms by the FIND software and data outside two standard deviations were discarded. 

The axial velocity (V z ) and tangential velocity (Vtheta) were combined to form the 
total velocity (V(otal) as shown in Figure 11. Flow velocities are seen to decrease toward the 
case wall. The boundary layer was distorted due to secondary flows and wake and comer 
vortices that form within the blade passage. Two-dimensional endwall flow angle was defined 
as the arc tangent of the tangential velocity divided by the axial velocity and is shown 
graphically in Figure 12. Periodicity is evident over 11.6 degrees (31 blades). 

Turbulence intensities were calculated with respect to the maximum downstream exit 
velocity (Vexit)- Figures 13 and 14 show turbulence intensity in both the tangential and axial 
directions. The tangential turbulence intensity is seen to be higher, possibly due to the steep 
gradient of the Vtheta mean-flow profile in the radial direction. 

The 0.9054 pressure ratio provided excellent seeding conditions and ATC vibrations 
were low. Data rates ranged from approximately 300 samples per second at the 3.34 mm 
depth to 150 samples per second at the 0.01 mm depth. Seeding material slowly accumulated 
inside the lower portion of the optical access hole and interfered with LDV data acquisition. 
Occasionally, the atomizer was secured and seed material allowed to disperse. Hole alignment 
and laser power (1.5 Watts) were periodically verified as these could drift due to temperature 
changes. Endwall measurement techniques did not provide pressure equalization across the 
optical access hole. During one exploratory run a minimum pressure ratio of 0.80 was 
achieved before the data rates deteriorated to unacceptable levels. 

A random mode comparison resulted in an average repeatability difference of 2.9%, 
0.6%, 6.4%, and 7.0% for velocity, flow angle, tangential turbulence intensity, and axial 
turbulence intensity respectively. Figures 15 through 23 graphically depict velocity, flow 
angle, and turbulence intensity random mode repeatability data for wake positions of +7, 0, 
and -8 degrees. A random and coincidence mode comparison resulted in an average 
maximum difference of 0.4%, 0.2%, 4.1%, and 6.3% for velocity, flow angle, tangential 
turbulence intensity, and axial turbulence intensity respectively. Tables FI and F2 of 
Appendix F contain LDV repeatability data. 



LDV Endwall Velocity 


Depth (mm) 



Figure 11. Two Component LDV Endwall Velocity 


LDV Endwall Flow Angle 



Figure 12. Two Component LDV Endwall Row Angle 


19 
























rbulence Intensity (%) 



Figure 13. Two Component LDV Endwall Turbulence, Vtheta 


V z Turbulence Intensity 

Referenced to Exit Velocity Depth (mm) 



Figure 14. Two Component LDV Endwall Turbulence, Vz 
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theta (m/s) 


Figure 15. LDV Endwall Flow Velocity Repeatability At +7 Degrees Wake Position 



Figure 16. LDV Endwall Flow Velocity Repeatability At 0 Degrees Wake Position 


LDV Endwall Flow Velocity 
Repeatability (-8° Wake) 



Depth (mm) 


Figure 17. LDV Endwall Flow Velocity Repeatability At -8 Degrees Wake Position 
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Figure 18. LDV Endwall Flow Angle Repeatability At +7 Degrees Wake Position 



Figure 19. LDV Endwall Flow Angle Repeatability At 0 Degrees Wake Position 



Figure 20. LDV Endwall Flow Angle Repeatability At -8 Degrees Wake Position 
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Figure 22 LDV Endwall Turbulence Repeatability At 0 Degrees Wake Position 


LDV Endwall Turbulence Intensity 
Repeatability (-8° Wake) 



Figure 23. LDV Endwall Turbulence Repeatability At -8 Degrees Wake Position 
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C. NUMERICAL COMPARISON 

Computational Fluid Dynamic (CFD) solutions were run with the same pressure ratios 
used in the experiments. Appendix G briefly describes the computational data reduction 
process and programs utilized. Graphical output of the numerical results was useful in 
obtaining a qualitative understanding of the flow characteristics. 

1. Blade Midspan Surface Pressures 

Reference 6 predicted hub, midspan and tip surface pressures at a pressure ratio of 
0.68. Figure 24 shows a comparison between numerical and experimental blade surface 
pressures at midspan for a pressure ratio of 0.6815. The suction surface curve suggests that 
the position of the blade passage throat was at 0.8 of axial chord and, in fact, the throat was 
located between static ports five and six. The minimum and maximum values at static ports 
five, six, and seven may be due to the boundary layer interaction at the blade's blunt trailing 
edge. Figures 25 through 28 show a comparison between numerical and experimental blade 
midspan surface pressures at pressure ratios of 0.5070, 0.6041, 0.8077 and 0.9054 
respectively. Comparison of the blade surface measurements and the computational results 
generally show excellent agreement. 
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0.5070 Pressure Ratio 
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Figure 25. P/Po vs. x/c for 0.5070 Pressure Ratio 
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Figure 26 P/Po vs. x/c for 0.6041 Pressure Ratio 
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0.8077 Pressure Ratio 



Figure 28. 




2. Mach Number and Flow Angle 

Two-dimensional Mach number and flow angle comparisons were made between 
numerical and experimental results for the surveys at each of the seven radial positions. 
Figures 29 through 32 compare Mach number, while Figures 33 through 36 compare flow 
angle at four selected radial positions. The remaining comparisons are included in Appendix 
H. All comparisons were based on a circumferential match between the experiment and 
computation, by shifting the latter to coincide with the measured profile. The specific 
matching was based on the maximum Mach number at the deepest radial position (3.34 mm). 
This matching was kept constant for the circumferential comparisons at all other spanwise 
locations. To depict flow periodicity and equivalently compare with experimental wake 
positions, the numerical solutions were repeated over one and one-half blade passages. 

Computed Mach numbers differed from the experimental values by an average 
12.7% between the surveys at 0.18 mm and 3.34 mm. The difference increased to 70.8% 
near the end wall, possibly due to a combination of LDV data velocity biasing and insufficient 
grid resolution. Velocity biasing was estimated to cause an approximate 11.0% velocity 
increase in the endwall region. The numerical solution at the endwall (0.01 mm) was based 
on the last k grid point in the radial direction. Computed flow angles surprisingly only varied 
from the experimental values by an average 3.2% (2.3 degrees). This was in contrast to the 
measurements in Reference 6 (Prat = 0.68) which showed poor comparison between 
measured and computed flow angles across the wake at 90 percent span. 

The question of radial spatial error with the probe volume needed to be determined. 
The probe volume length and diameter (Figure 37) was calculated as 1.56 mm and 0.11 mm 
respectively. The probe volume had approximately 30 fringes across its minor axis and after 
data processing an average of 13 fringe crossings constituted a valid Doppler burst. Since the 
IFA-750 digital burst correlator automatically centered the valid signal within the probe 
volume, an effective probe volume with length 1.39 mm and diameter 0.04 mm was obtained. 
Figure 38 depicts, to scale, the radial survey resolution that results from endwall depth and 
effective probe volume dimensions. The decrease in relative spatial resolution close to the 
endwall also contributed to the lack of agreement between the Mach number profiles (Figure 
32). 
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Figure 30. Mach Number Comparison At 0.89 mm Depth 
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Figure 31. Mach Number Comparison At 0.18 mm Depth 


Mach Number Comparison 
0.01 mm Depth (99.9% Span) 
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Figure 32. Mach Number Comparison At 0.01 mm Depth 
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Wake Flow Angle Comparison 
3.34 mm Depth (80.9% Span) 
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Figure 34. Wake Flow Angle Comparison At 0.89 mm Depth 
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Figure 35. Wake Flow Angle Comparison At 0.18 mm Depth 
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Effective Probe Volume with 13 Fringes 












D. NUMERICAL RESULTS 

I. Convergence History 

Solutions at pressure ratios of 0.5070, 0.6041 and 0.6815 fully converged after 
approximately 1,200 iterations. Figure 39 shows density residuals at 0.6815 pressure ratio 
after 10,200 iterations. Solutions for pressure ratios of 0.8077 and 0.9054 were initially 
obtained after 1,200 iterations, but convergence did not occur (i.e. three orders of magnitude 
reduction of the residuals) until approximately 3,000 iterations. Figures 11 through 15 of 
Appendix I contain density-residual convergence history for all examined pressure ratios. 


RVC3D Residuals [0.6815 Pressure Ratio] 

1.00E-04 

w 1.00E-05 
| 1.00E-06 
•3 1.00E-07 


1.00E-08 

1 2000 4000 6000 8000 10000 

Number of Iterations 


Figure 39. 0.6815 Pressure Ratio Convergence History (10,200 Iterations) 

2. Turbulence Models 

At a pressure ratio of 0.8077, the code was run with the following turbulence models; 
RVC3D version 206 Cebeci-Smith 

RVC3D version 206 Baldwin-Lomax 
RVC3D version 208 Baldwin-Lomax 

so as to assess the effect of different turbulence models on the solution. In all three cases, the 
code produced similar solutions. 
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3. Mach Number Contours 

Figures 40 and 41 display exit plane Mach number contours at pressure ratios of 
0.9054 and 0.5070 respectively. The Mach number contours between 25% and 75% span 
provide an indication of where the flow characteristics were well established. Boundary layers 
appear, as expected, to be thin near the endwall region and thicker near the hub. Spanwise 
wake curvature, due to secondary flows, is also apparent in the two figures with Figure 40 
showing the greater curvature. 

Figure 42 displays midspan Mach number contours for the blade passage at 0.5070 
pressure ratio. The flow was strongly accelerated along the suction side. The acceleration of 
the flow to high local Mach numbers resulted in shocks at the throat and trailing edge. The 
wake was a merging of the pressure and suction side boundary layers. The trailing-edge 
shock extended across the wake from the adjacent (upper) blade, out to the exit plane of the 
grid. The shock interaction with the wake caused the wake to narrow down, resulting in a 
wake which was repeatedly diffused and then coalesced. The shock strength decreased away 
from the trailing edge as a result of this interaction with the wake, which could be both 
physical and computational in nature. Because of the relatively good agreement between 
experiment and computation of the blade surface pressure field, it was felt that the solution 
was realistic. Appendix J includes program listings for a graphical display of multiple grids 
and corresponding numerical solutions for use with PLOT3D. 




Figure 40. Exit Plane Mach Contours (0.9054 Pressure Ratio) 
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Figure 41. Exit Plane Mach Contours (0.5070 Pressure Ratio) 
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V CONCLUSIONS AND RECOMMENDATIONS 


An experimental and computational investigation was conducted of the flow through 
an annular turbine cascade. Computational pressure ratios (Prat), defined as the downstream 
hub-static pressure (Phub) divided by the upstream stagnation pressure (Po). were chosen to 
coincide with experimental test conditions. Downstream computational and experimental 
measurement locations were chosen to coincide to provide a proper comparison. 

An existing annular turbine cascade was modified to provide a capability to measure 
blade midspan surface pressures. Experimental measurements of blade midspan surface 
pressures were obtained and favorably compared with numerical predictions over a range of 
subsonic and transonic conditions. Measurements were repeatable with a maximum 
repeatability error of less than 0.87%. Numerical solutions converged on pressure within 
approximately 1,200 iterations. 

The minimum and maximum pressure values on the suction surface close to the 
trailing edge should be investigated further to confirm the influence of the relatively blunt 
trailing edge. An inlet boundary layer survey, requiring a rig modification, should be 
performed to supplement data in Reference 6 and refine the inlet boundary conditions to the 
computational model. 

To supplement blade midspan surface pressures, blade pressure taps could be installed 
at other spanwise locations to include the hub, tip, and possibly end wall locations. A tip gap 
could be created and a computational tip gap model introduced to investigate the endwall 
region. Blade tips, and the adjacent endwall, could be instrumented for pressure readings and 
a numerical and experimental comparison conducted to obtain insight into tip leakage flows. 
In addition, an LDV optical-access window could be installed over the tip region to measure 
tip-gap flow characteristics. 

Successful measurements of two-dimensional velocity, flow angle, and turbulence 
intensity through a 1.0922 mm diameter casing-access hole were demonstrated. All LDV 
data were acquired one-half axial chord downstream over seven radial positions and 17 
separate peripheral displacements. Experimental flow velocities were increased from those of 
Reference 6 mainly due to the availability and use of a digital burst correlator in place of a 
signal processor. Although all LDV data reported herein were acquired at 0.9054 pressure 
ratio, measurements were attempted at 0.80 pressure ratio and some data were obtained; 
however, low data rate and endwall hole vibrations prevented a complete survey. 
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The comparison of computational and experimental Mach number and flow angle 
yielded an average 12.7% difference in Mach number and 3.2% (2.3 degrees) difference in 
flow angle. All experimental measurements were repeatable to an average uncertainty 
ranging from 0.2% to 7.0% for velocity, flow angle, and turbulence intensity in both random 
and coincidence data processor modes. All numerical and experimental comparisons were 
based on a circumferential match of maximum Mach number at the 80.9% (3.34 mm) span 
location. An investigation of other circumferential matching options should be conducted to 
explore the effect on the degree of agreement (i. e. minimum Mach number circumferential 
matching). 

LDV alignment and radial survey procedures were devised and repeated. The probe 
volume passed cleanly through the optical access hole and radial surveys were conducted to 
provide complete coverage for data recording. Probe volume dimensions were further 
refined to aid in estimates of radial-position accuracy. Introduction of a second fiber-optics 
probe for three component measurements would decrease the effective probe volume and 
increase radial position accuracy. The ability to resolve endwall flow characteristics and 
eventually flow tip losses appears to be promising. 

The present hardware did not provide pressure equalization across the optical access 
hole. A rig modification is required to equalize the pressure at the measurement location and 
eliminate any influence the opening might have on the flow characteristics. Also, the coarse 
peripheral (wake) positioning mechanism needs modification to allow more precise 
adjustment. Experimental peripheral positioning uncertainties may contribute to differences 
in the compared data. 

Future efforts with computational fluid dynamics should include modeling the inlet 
bellmouth and center body, and conducting a turbulence model sensitivity analysis in the 
endwall region. Increased grid resolution in the endwall region would be beneficial. The 
experimental 99.9% span location coincided with the final k grid point and may account for 
the 70.8% difference between the experimental and computed Mach numbers. Finally, a 
detailed investigation into transonic and sonic flow conditions could be conducted to examine 
the blade-passage and trailing-edge shock locations and subsequent interactions, both 
computationally and experimentally. 
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APPENDIX A. PRESSURE DATA ACQUISITION 

Pressure data acquisition was accomplished using a Hewlett-Packard 9000 computer 
system with the following program: 

10 1FILE NAME: BILLSP 

20 IDISK LABEL: '"/AUSTIN" 

30 1LAST MODIFIED 10/27/94 BY LT BILL DONOVAN 

40 !THIS PROGRAM RECORDS AND REDUCES MEASURED PRESSURES FROM A 

5 0 1SCANIV ALVE CONNECTED TO 11 STATIC PRESSURE PORTS ON ATC 

60 ! VARIABLES / PORT CONNECTIONS... 

70 ! V=DESIRED S/V 

80 1 A=PRESENT S/V PORT 

90 ! C=CHANNEL 

100 ! Voltage ()=DVM VOLTAGE READING FOR EACH PORT 

110 ! S=SCANNER NUMBER 

120 ! Pinhg=BAROMETRIC PRESSURE IN INHG 

130 ! Pamb=BAROMETRIC PRESSURE (psi) 

140 ! AMBIENT PRESSURE...# 1 

150 ! CALIBRATION PRESSURE...#2 

160 ! REFERENCE UPSTREAM TOTAL PRESSURE...PORT #4 

170 ! INNER HUB STATIC PRESSURE...PORT #5 

180 ! UPSTREAM STATIC PRESSURE...PORT #6 

190 ! ATC PASSAGE PORTS #1-#11 ...S/V PORTS #9-#19 

200 ! Prat= INNER HUB STATIC/ REFERENCE UPSTREAM TOTAL 

210 ! Tstag=UPSTREAM STAGNATION TEMPERATURE IN DEGREES F 

220 OPTION BASE 1 

230 CLEAR SCREEN 

240 PRINTER IS CRT 

250 DISP "PLEASE WAIT WHILE RESETTING SCANIVALVE TO PORT #1" 

260 PRINT 
270 V=5 

280 A=1 

290 GOSUB Read 
300 DIM Y(20) 

310 DIM Voltage(50) 

320 DIM P(50) 

330 INPUT "ENTER MONTH, DAY, YEAR (I.E.02, 04, 94)", Y(3), Y(5), Y(7) 

340 INPUT "ENTER RUN #: ", Y(9) 

350 PRINTER IS 711 

360 PRINT USING ”K, DD, " ” / " ", DD, " " / " ", DD", "DATE OF RUN: ", Y(3), Y(5), Y(7) 
370 PRINT "DATA RUN ", Y(9) 

380 PRINT USING 7,5X, K, DDD,/, 2X, K, 15X, K, 17X, K", "SCANIVALVE# ", V, 
"Pott", "Volts", "Psi” 

390 PRINTER IS CRT 

400 PRINT " ZERO DVM ON SCANIVALVE #5, PORT #1..." 

410 PRINT 

420 PRINT " SET 20 inhg FOR SCANIVALVE CALIBRATION..." 

430 PRINT 

440 PRINT " ADJUST SPAN FOR DVM = 0.009823 VOLTS ON PORT #2..." 

450 PRINT 

460 PRINT "ENSURE DESIRED PRESSURE RATIO IS SET..." 

470 INPUT "ENTER BAROMETRIC PRESSURE IN INHG", Pinhg 

4 80 INPUT "ENTER HIGHEST SCANIVALVE PORT # NEEDED ”, Nports 

490 INPUT "ENTER UPSTREAM STAGNATION TEMP. IN DEG. F ", Tstag 





500 Pamb=Pinhg*(0.4911541) 

510 PRINT "**** PRESS CONTINUE WHEN READY TO TAKE DATA***" 

520 PAUSE 

530 ,**„.******************************************„»***,„* 

540 !RECORD THE PRESSURE DATA FOR THIS RUN 

550 OUTPUT 722; "F1R1T1M0Z1" ! SETS UP HP3456A DVM 

560 V=5 ISCANIVALVE #5 DESIGNATED 

570 S=1 !HP3495A SCAN #1 DESIGNATED 

580 FOR A=1 TO Nports 

590 GOSUB Read 

600 WAIT 2.0 

610 OUTPUT 701; "C" !CLEARS SCANNER #1 

620 OUTPUT 701 USING "DDD"; V+9 [SCAN SET TO READ S/V 

630 TRIGGER 722 

640 ENTER 722; Voltage(A) 

650 P(A)=(Voltage(A) - Voltage(l))*1000+Pamb 

660 PRINTER IS CRT 

670 PRINT USING ”2X,DDD,7X,7DD.DDDDDDD.8X,7DD.DDDD"; A, Voltage(A), P(A) 
680 PRINTER IS 711 

690 PRINT USING "2X,DDD,7X,7DD.DDDDDDD,8X,7DD.DDDD"; A, Voltage(A), P(A) 
700 NEXT A 

710 Prat=P(5)/P(4) 

720 PRINT "Prat = ", Prat 

730 PRINT "UPSTREAM STAGNATION TEMPERATURE WAS (DEG F)", Tstag 

740 PRINT "INPUT BAROMETRIC PRESSURE IN (inhg) WAS ”, Pinhg 

750 !**********SUBROUTINE TO POSITION AND READ S/V PORTS********* 

760 Read: OUTPUT 707 USING "#,K";V 
770 PO=SPOOL(707) 

780 L=BINAND(P0,15) 

790 T1=ROTATE(PO,4) 

800 M1=BINAND(T1,7) 

810 P(A)=10*M1+L 

820 CLEAR 707 

830 IF P(A)=A THEN Finish 

840 OUTPUT 701; "C" 

850 OUTPUT 701 USING "DDD";V-1 

860 OUTPUT 701; "C" 

870 WAIT 0.1 

880 GOTO Read 

890 Finish: RETURN 
900 END 
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The following table relates each Scanivalve port to its respective pressure 
measurement. 


Scanivalve Port Number 

Pressure Measured 

1 

Ambient Pressure 

2 

Calibration Pressure (set at 20 inches hg) 

3 

Not Used 

4 

Upstream Total Pressure (Po) 

5 

Downstream Hub Static Pressure (Phub) 

6 

Upstream Static Pressure (Pstatic) 

7 

Not Used 

8 

Not Used 

9 

Blade Static Port #1 (at leading edge) 

10 

Blade Static Port #2 (suction side) 

11 

Blade Static Port #3 (suction side) 

12 

Blade Static Port #4 (suction side) 

13 

Blade Static Port #5 (suction side) 

14 

Blade Static Port #6 (suction side) 

15 

Blade Static Port #7 (suction side) 

16 

Blade Static Port #8 (at leading edge) 

17 

Blade Static Port #9 (pressure side) 

18 

Blade Static Port #10 (pressure side) 

19 

Blade Static Port #11 (pressure side) 


Table Al. Pressure Data Acquisition Connections 
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APPENDIX B. GRID GENERATION INPUT FILE 


Grid generation was completed on Silicon Graphics Indigo II workstations. The 
following TCGRID namelist file was utilized to obtain the final grid. 


&naml im=150 jm=31 km=65 itl=20 icap=18 k2d=3 merid=0 &end 
&nam2 nle=16 nte=10 dsle=.018 dste=.003 dshub=.0004 dstip=.0004 
dswte=.001 dswex=.060 dsthr=l, dsmin=.0004 dsmax=.025 
dsra=.45 rcom=.098 &end 

&nam3 iterm=100 idbg=0 0 0 0 0 0 0 0 aabb=.5 ccdd=.45 &end 
&nam4 zbc=-1.5 -1.5 0.5 -1.5 -1.5 .5 

rbc= 3.895 3.895 3.895 4.585 4.585 4.585 &end 
TRANSONIC TURBINE' 

2 2 

-.975 0.0 
3.895 3.895 
-.975 0.0 
4.585 4.585 
2 51 31 

-7.8999996E-03 -9.920000IE-03 -1.2000000E-02 -1.4080000E-02 
-1.6100001E-02 -1.7999999E-02 -1.971000IE-02 -2.1190001E-02 
-2.2390001E-02 -2.3280000E-02 -8.1040002E-02 -0.1503800 
-0.2308900 -0.3219400 -0.4228500 -0.5328500 -0.6511000 
-0.7893000 -0.9275000 -0.9410600 -0.9527700 -0.9622700 
-0.9692700 -0.9735600 -0.9750000 -0.9735600 -0.9692700 
-0.9622700 -0.9527700 -0.9410700 -0.9275000 -0.9124900 
-0.8250000 -0.7200000 -0.6000000 -0.4560000 -0.3440000 
-0.2730000 -0.2200000 -0.1820000 -0.1160000 -5.9000000E-02 
-1.7000001 E-02 O.OOOOOOOE+OO -1.8220000E-04 -7.2359998E-04 
-1.6100000E-03 -2.8100000E-03 -4.2900001E-03 -6.000000IE-03 
-7.8999996E-03 -0.2349885 -0.2351271 -0.2351733 -0.2351271 
-0.2349885 -0.2347600 -0.2344519 -0.2340719 -0.2336329 
-0.2331451 -0.1980565 -0.1643980 -0.1324210 -0.1023671 
-7.4464694E-02 -4.8926830E-02 -2.5946086E-02 -4.6726577E-03 
1.6618744E-02 1.9057767E-02 2.2061616E-02 2.5545571E-02 
2.9399229E-02 3.3504494E-02 3.7740692E-02 4.1976891E-02 
4.6082158E-02 4.9935814E-02 5.3417202E-02 5.6423619E-02 
5.8862645E-02 6.0659818E-02 6.9319643E-02 7.1887039E-02 
6.7522466E-02 5.1347882E-02 2.5673941E-02 O.OOOOOOOE+OO 
-2.5673941E-02 -5.1347882E-02 -0.1026958 -0.1540437 
-0.2053915 -0.2320924 -0.2326265 -0.2331451 -0.2336329 
-0.2340719 -0.2344519 -0.2347600 -0.2349885 3.895000 

3.895000 3.895000 3.895000 3.895000 3.895000 3.895000 

3.895000 3.895000 3.895000 3.895000 3.895000 3.895000 

3.895000 3.895000 3.895000 3.895000 3.895000 3.895000 

3.895000 3.895000 3.895000 3.895000 3.895000 3.895000 

3.895000 3.895000 3.895000 3.895000 3.895000 3.895000 

3.895000 3.895000 3.895000 3.895000 3.895000 3.895000 

3.895000 3.895000 3.895000 3.895000 3.895000 3.895000 

3.895000 3.895000 3.895000 3.895000 3.895000 3.895000 

3.895000 3.895000 -7.8999996E-03 -9.920000IE-03 
-1.2000000E-02 -1.4080000E-02 -1.6100001E-02 -1.7999999E-02 
-1.9710001 E-02 -2.1190001 E-02 -2.2390001 E-02 -2.3280000E-02 
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-8.1040002E-02 -0.1503800 -0.2308900 -0.3219400 -0.4228500 
-0.5328500 -0.6511000 -0.7893000 -0.9275000 -0.9410600 
-0.9527700 -0.9622700 -0.9692700 -0.9735600 -0.9750000 
-0.9735600 -0.9692700 -0.9622700 -0.9527700 -0.9410700 
-0.9275000 -0.9124900 -0.8250000 -0.7200000 -0.6000000 
-0.4560000 -0.3440000 -0.2730000 -0.2200000 -0.1820000 
-0.1160000 -5.9000000E-02 -1.700000IE-02 O.OOOOOOOE+OO 
-1.8220000E-04 -7.2359998E-04 -1.6100000E-03 -2.8100000E-03 
-4.290000IE-03 -6.0000001 E-03 -7.8999996E-03 -0.1996249 
-0.1997426 -0.1997819 -0.1997426 -0.1996249 -0.1994308 
-0.1991690 -0.1988462 -0.1984733 -0.1980589 -0.1682508 
-0.1396576 -0.1124929 -8.6961828E-02 -6.3258447E-02 
-4.1563794E-02 -2.2041440E-02 -3.9694658E-03 1.4117776E-02 
1.6189748E-02 1.8741548E-02 2.1701200E-02 2.4974918E-02 
2.8462378E-02 3.2061070E-02 3.5659760E-02 3.9147221E-02 
4.2420939E-02 4.5378406E-02 4.7932386E-02 5.0004359E-02 
5.1531076E-02 5.8887679E-02 6.1068702E-02 5.7360962E-02 
4.3620501E-02 2.1810250E-02 0.0000000E+00 -2.1810250E-02 
-4.3620501E-02 -8.7241001E-02 -0.1308615 -0.1744820 
-0.1971647 -0.1976183 -0.1980589 -0.1984733 -0.1988462 
-0.1991690 -0.1994308 -0.1996249 4.585000 4.585000 4.585000 

4.585000 4.585000 4.585000 4.585000 4.585000 4.585000 

4.585000 4.585000 4.585000 4.585000 4.585000 4.585000 

4.585000 4.585000 4.585000 4.585000 4.585000 4.585000 

4.585000 4.585000 4.585000 4.585000 4.585000 4.585000 

4.585000 4.585000 4.585000 4.585000 4.585000 4.585000 

4.585000 4.585000 4.585000 4.585000 4.585000 4.585000 

4.585000 4.585000 4.585000 4.585000 4.585000 4.585000 

4.585000 4.585000 4.585000 4.585000 4.585000 4.585000 



Figure Bl. Blade Geometry From Ref. [6] 
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APPENDIX C. RVC3D SAMPLE INPUT AND EXECUTION 


The following input namelist example was for a 0.6815 pressure ratio with subsonic 
exit boundary conditions and Cebeci-Smith turbulence modeling. 


TRANSONIC TURBINE Annular Cascade' 

&nll im=150 jm=31 km=65 itl=20 iil=67 &end 

&nl2 cfl=5.5 avisc 1=0.0 avisc2=0.0 avisc4=0.50 ivdt=l nstg=4 itmax=1200 
irs=l epi=0.50 epj=0.60 epk=0.60 &end 
&nl3 ibcin=3 ibcex=3 isymt=0 ires=10 icmt=50 
iresti=0 iresto=l ibcpw=0 iqin=0 &end 
&nl4 emxx=0.13 emty=0.0 emrz=0.0 expt=0.0 prat=0.6815 ga=1.4 
om=0.000000 igeom=l alex=-67.0 &end 
&nl5 ilt=3 tw=1.00 renr=6.651e6 pmr=.7 prtr=.9 vispwr=.666666 

srtip=0.0 cmutm=14. jedge=31 kedge=20 iltin=2 dblh=0.0048 dblt=0.0178 &end 
&nl6 iol=l io2=165 oar=0. ixjb=0 njo=l nko=3 
jo=l ko=5 11 16 &end 


Runs on the NPS Cray computer took approximately 12 hours for a 1,200 iteration 
convergence and were accomplished with the command: 

>qsub myjob.nqs 

where the file myjob.nqs was as follows: 
myjob.nqs 

#QSUB -q prem -IT 50000 -1M 9Mw 
cd /dl/whdonova 

./exec.rvc3d < rvc3d.in > rvc3d.out 


Runs on the NASA Ames Research Center Cray computer took approximately 1 hour 
for a 1,200 iteration convergence and were accomplished with the command: 

>qsub -lm lOMw -It 14400 rvc3d.com 
where the file rvc3d.com was as follows: 
rvc3d.com 

cd /m/vn/v3/wdonovan 
,/exec.rvc3d < rvc3d.in > rvc3d.out 
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APPENDIX D. BLADE MIDSPAN SURFACE PRESSURE DATA 
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Table Dl. 0.6815 Pressure Ratio Series #1 
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Table D2. 0.6815 Pressure Ratio Series #2 
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Table D3. 0.5070 Pressure Ratio Series 
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Table D4. 0.6041 Pressure Ratio Series 
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Table D5. 0.8077 Pressure Ratio Series 
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Table D6. 0.9054 Pressure Ratio Series 
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APPENDIX E. LDV DATA 


Wake Position 

OepUt 

v-theta 

V-* 

V tout 

Mach 

t low Angle 

v-theu 


1 stag (“P-) 

95.0 

(degrees) 

(mm) 

(m/s) 


(m/s) 

.Number 

(degrees) 

turb.int. 

turb. int 

P static (inhg) 

34.63 

-8 






/53 

U6 

5.19 

P slag (Inhg) 

3.i3 











0.9056 

-8 









M exit 

"6574 












•a 









a exit (m/s) 

346.84 

-8 







5.41 

3.98 

a stag (m/s) 

351.79 









3.98 

V exit (m/s) 

131.46 












Wake Position 

D cp [ ii 




Mach 




T stag ( J b‘) 

-srr 

(degrees) 

(mm) 


(m/s) 


Number 

(degrees) 

turb.int 

turb. int 

P static (inhg) 









756 


P slag (inhg) 


- 1 









Prat 

0.9056 








7.83 


M exit 

0579 











299.4 



96.6 




/6.6 

6.96 


a exit (m/s) 

346.84 

■' 






76.8 

3.49 

4.18 

a stag (m/s) 


: 






/6.2 

5.09 



131.46 












Wake Position 

Depth 

V-theta 

V-z 

V total 

Mach 

How Angle 

V-theta 


1 stag TO 

95.0 

(degrees) 

(mm) 

(m/s) 

(m/s) 

(m/s) 

Number 

(degrees) 

lurb. int. 

turb. int. 

P static (Inhg) 

30.03 

-6 









P stag (inhg) 


-6 


87.6 




76.0 

7.56 

5.66 

Prat 

0.9056 

■* 









Mextt 

0579 

-6 









Texit(V) 

299.4 

•6 






76.6 

6.40 


a exit (m/s) 

346.44 

•* 

1./8 

99.7 

23.5 

102.4 

0.291 

76.7 

526 

4.29 

a stag (m/s) 


6 






16.1 



V exit (m/s) 

131.46 












Wake Position 

iiepth 

V -theta 


V total 

Mach 

How Angle 

V-theta 

V-z 

1 stag TO 

96.0 

(degrees) 

(mm) 

(m/s) 

(m/s) 

(m/s) 

Number 

(degrees) 

turb. int. 

turb. int. 

P static (Inhg) 

30.03 

- 5 






15.1 


652 

P stag (inhg) 

3.i3 


0.06 





75.6 



Prat 

0.9056 








7.29 

5.90 


0579 

5 






/6.1 



T exit ( D K) 

299.9 

- s 



23.6 






a exit (m/s) 








163 

5.11 

4./5 

a stag (m/s) 


‘ 5 

3.34 

102.0 

246 

104 9 

0.298 

I6A 

5.9/ 

4.13 




Table El. LDV Data (Wake Positions of -8, -7, -6, and -5 Degrees) 
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Table E2. LDV Data (Wake Positions of -4, -3, -2, and -1 Degrees) 
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Wake Position 

Depth 

V-theta 



Mach 

How Angle 



T stag ("k, 

96.0 

(degrees) 





Number 




P static (inhg) 




90.2 

24./ 



M.I 



P stag (inhg) 



0.06 






7.68 

4.89 

Prat 

0.9056 


0.18 







4.68 

M ezit 

6.37$ 










1 exit ("K) 

299.9 

0 

0.89 





/6.5 

6.10 

4.13 

a exit (m/s) 




108.6 




/6.5 

3.90 


> stag (m/s) 


U 


106.8 




















Make Position 

Depth 

V-thela 



Mach 

How Angle 

V-theta 


1 stag (°K) 

100.0 

(degrees) 


(m/s) 

(m/s) 


Number 

(degrees) 

turb. ml. 

turb. int. 

P static (inhg) 

36.63 








8.18 


P stag (inhg) 



0.06 







4.63 

Prat 

0.9056 







/5.6 

8.03 

4.59 

M exit 

0J79 


6.4i 







4.11 


342.1 







/6.4 


3-77 

a exit (m/s) 



l./S 






3.96 

3.68 

a stag (m/s) 


1 




110.3 

0.312 

/6.0 


3.6/ 

V exit (m/s) 

132.0a 












Wake Position 



V-z 


Mach 

How Angle 



1 stag TO 

100.0 

(degrees) 


(m/s) 


(m/s) 

Number 

(degrees) 

turb. int. 

turb. ml. 

P static (inhg) 

36.63 










P stag (inhg) 


2 











1 







8.18 

4.45 

M exit 

0-379 

2 









1 exit (°K) 

302.1 




23.6 






a exit (m/s) 

346.46 

2 









a stag (m/s) 

353.37 









3.66 

V exit (m/s) 

131.05 












Make Position 

Depth 

V-theta 

V-z 


Mach 

I*low Angle 

V-theta 

V-z 

T stag TO 

9X0 

(degrees) 

(mm) 

(m/s) 

(m/s) 

(m/s) 

Number 

(degrees) 

turb. int. 

turb. ml. 

P static (inhg) 

30.03 

3 





0.261 



4.74 

P stag (inhg) 



TT55 - 







4.64 

Prat 

0.9056 

3 

0.18 



94.7 


76.0 


4.01 

Mexit 

0-379 







7S!5 

7.66 

3.94 

T exit (°K) 

301.0 










a exit (m/s) 









5.28 

3.18 

a stag (m/s) 





25.6 






V exit (m/s) 



Table E3. LDV Data (Wake Positions of 0, 1,2, and 3 Degrees) 
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Table E4. LDV Data (Wake Positions of 4,5, 6, and 7 Degrees) 
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Table E5. LDV Data (Wake Position of 8 Degrees) 
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APPENDIX F. LDV REPEATABILITY DATA 


Wake Position 

Dt P lh 

V .theu 

v-* 

V total 

Mach 

1- low Angle 


■ V-z 

'I'stag ("10 

90.0 

(degrees) 


(m/s) 

(m/s) 

(m/s) 

Number 

(degrees) 

turb. ink 

turb. ml. 

1 

29.94 

Coincidence -8 







8./4 

4.61 

P slag (mhg) 

3.13 

Coincidence -8 

0.05 








Prat 

0,9055 








8.66 

4.54 

M exit 

0J79 

Coincidence -6 




98.9 




4.48 


299.9 

Coincidence -it 

0.89 





nA 

7.58 

4.27 

a exit (m/s) 

J47.15 

Coincidence -S 

1.78 



104.5 





a stag (m/s) 

551.16 

Coincidence -8 





0.304 




V exit (m/s) 













Wake Position' 


V-theta 


V tout 

Mach 

(■low Angle' 

V-theta 


t stag HO 

96.0 

(degrees) 

(mm) 

(m/s) 

(m/s) 

(m/s) 

Number 

(degrees) 

turb. ink 

turb. ink 

P static (inhg) 

30.00 

Random -8 





0.263 

/6.0 


4-'4 

P stag (inhg) 

3.13 

Random -8 


91.4 






4.67 


6.4655 

Random -8 









M exit 


Random -8 








4.44 



Random -8 









a exit (m/s) 


Random -8 









a stag (m/s) 

352.10 

Random -8 







5.16 


V exit (m/s) 

131.58 












Wake Position 

Depth 


V-z 


Mach 

Mo" Angle 


v-z 

1 stag (°P) 

556 

(degrees) 

(mm) 

(m/s) 


(m/s) 


(degrees) 

turb. ink 

turb. ink 

P static (inbg) 

30.00 

Coincidence 0 









P stag (inhg) 


Coincidence 0 

0.06 








Prat 

0.9055 

Coincidence 0 





0276 



5.14 

M exit 

0379 

Coincidence U 








4-59 



Coincidence 0 






16.6 

6.58 


a exit (m/s) 


Coincidence 0 






16.6 



a stag (m/s) 


Coincidence 6 









V exit (m/s) 













Wake Position 

Depth 




Mach 

(•low Angie 



T stag HO 

87 0 

(degrees) 

(mm) 

(m/s) 

(m/s) 

(m/s) 

Number 

(degrees) 

turb. int. 

turb. ml. 

Pa h 


Random 0 









P stag (inhg) 

3.13 

Random 0 

0.06 








Prat 

0.9055 

Random 0 

O.iS 







4.64 

M exit 


Random 0 






76.1 

6 .87 


1 exit (°K) 

295.1 

Random t) 






7S3 

5./9 

3.82 

a exit (m/s) 

54433 

Random 0 






76.5 

4.32 


a stag (m/s) 

34934 

Random (J 



26.6 

110J 

UJI6 




V exit (m/s) 



Table FI. LDV Repeatability Data (Wake Positions of -8 and 0 Degrees) 
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I Random 7 | 134 | 103.0 | K.3 | IK. I | a.m ["" Ki | 4.18 | 4.50 | V e«ii fa/si | | 

Table F2. LDV Repeatability Data (Wake Position of 7 Degrees) 
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APPENDIX G. COMPUTATIONAL DATA REDUCTION 


Residual (fort.4) and P/Po (fort.7) files were produced from the solution (fort.3) files 
using the FORTRAN program pxy.f. 


c pxy.f reads rvc3d files & writes ascii files for gnuplot 
c unit 1 = input xyz file 
c unit 3 = input q file 

c unit 7 - output blade pressures on 5 k-planes 
c unit 4 = output residual history 

parameter(ni=150,nj=31 ,nk=65) 

real x(ni,nj,nk),y(ni,nj,nk),z(ni,nj,nk) 
real qq(5,ni,nj,nk),resd(5000,5) 
real pk(5),xk(5) 
dimension kk(5) 

c k-values are hard-wired below (hub, 25%, 50%, 75%, tip) 
data kk/2,27,33,39,64/ 

c read grid coordinates 
read(l,*)im,jm,km 

read(l,*)(((x(ij,k),i=l,im),j=ljm),k=l,km), 

1 (((y(i.j,k).>= 1 .im),j= 1 jm),k= 1 ,km), 

2 (((z(ij,k),i=l.im),j=l,jm),k=l,km) 

c read restart file 


icheck=iabs(im-imax)+iabs(jm-jmax)+iabs(km-kmax) 

if(icheck.ne.O)then 

write(6,610)im,jm,km,imax,jmax,kmax 

stop 

endif 

read(3,*)((((qq(l,i,j,k),i=l ,im) j=l jm),k=l,km),1=1,5) 

c additional residual data 

read(3,*)itl,iil,phdeg,ga,om,nres,dum,dum,dum,dum 
read(3,*)((resd(nr,l).ru^l ,nres),l=l ,5) 

c ps/pO output to unit 3 

itr=im+l-itl 

ggm=ga*(ga-l.) 

j=l 
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c normalize x by chord 
do 71=1,5 
k=kk(l) 

xmin=x(im/2,j,k) 
xmax=x(itl, j,k) 
do 5 i=itl,itr 
xmin=min(xmin,x(i,j Jc)) 

5 xmax=max(xmax,x(i,j,k)) 
chord=xmax-xmin 
do 7 i=itl,itr 

7 x(i,j,k)=(x(ij,k)-xmin)/chord 

write(7,305)(kk0)4= 1,5) 
do 20 i=itl,itr 
do 101=1,5 
k=kk(l) 

pk(l)=ggm*(qq(5,i,j,k)-.5*(qq(2,ij,k)**2+qq(3,i,j,k)**2 

l+qq(4,i,jjc)**2)/qq(l,i,jjc)) 

10 xk(l)=x(ij,k) 

20 write(7,300)i,(xk(l),pk(l),l=l ,5) 



write(4,310) l,(resd( 1,1)4= 1,5) 
do 30 j=2,nres 
it=10*(j-l) 

30 write(4,310)it,(resd(j,l),l=l,5) 

300 format(i5,10fl0.5) 

305 formate k=',5(17x,i3)) 

310 format(i5,5(lx,el0.3)) 

610 formatC ***** warning ****** j t 

1 ’ im, jm, km, read from input',3i5,' do not match',/, 

2 ' im, jm, km, read from restart file',3i5) 
stop 

end 
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The plane.f program from Reference 6 was modified to compare two-dimensional 
experimental and computational Mach number and wake flow angle. Mach number, based 
on stagnation conditions, was redefined with respect to tangential and axial velocity 
components. 


Program plane.f 

c Modified by Lt. William Donovan on 08 May 95 for Master’s Thesis 
c plane.f reads rvc3d files & writes ascii files for plotting 
c exit plane mach number and flow angles 
c (2D comparison with experimental results) 
c unit 1 = input xyz file 
c unit 3 = input q file 
c unit 4 = output residual history 

c unit 7 = output (k #1) counter, theta, pt/pO,Mach, Angle 

c unit 8 = output (k #2) counter, theta, pt/pO,Mach, Angle 

c unit 9 = output (k #3) counter, theta, pt/pO,Mach, Angle 

c unit 10 = output (k #4) counter, theta, pt/pO,Mach, Angle 

c unit 11 = output (k #5) counter, theta, pt/pO,Mach, Angle 


parameter(ni= 150,nj=31 ,nk=65) 
integer um 

real x(ni,nj,nk),y(ni,nj,nk),z(ni,nj,nk), gama 
real qq(5,ni,nj,nk),resd(5000,5),ql ,q2,q3,q4,q5,m(5) 
dimension kk(5),ang(5),theta(5),n0(5),vtnd(5),vmd(5), 
1 vznd(5),ptp0(S),deg(5) 


c k-values are hard-wired below (80.9%, 89.8%, 94.9%, 97.6%, 99.0%) 
c (use span.f to obtain %-k equivalency) 

C (k values of 61 and 64 obtained from a second program run) 
c 

data kk/41,45,49,53,57/ 
c read grid coordinates 

c ************************************************* 
read(l ,*)im,jm ,km 

read(l ,*)(((x(i j,k),i= 1 ,im),j= 1 jm),k= 1 ,km), 

1 (((y(i j,k),i= 1 ,im)j= 1 ,jm),k= 1 ,km), 

2 (((z(i,j,k),i=l ,im),j=l jm),k= 1 ,km) 


read restart file 


icheck=iabs(im-imax)+iabs(jm-jmax)+iabs(km-kmax) 

if(icheck.ne.O)then 

write(6,610)im,jm,km,imax,jmax,kmax 

stop 

endif 
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read(3,*)((((qq(l I i,j,k),i=l,im),j=l,jm),k=l,km),l=l,5) 


c additional residual data 


read(3,*)itl,iil,phdeg,gama,om,nres,dum,dum,dum,dum 
read(3,*)((resd(nr,l),nr= 1 ,nres),l=l ,5) 



k=kk(3) 

i=l 

do 10 j=jmax,l,-l 
um=um+l 
do 151=14 
k=kk(l) 
c 

c Rename desired q values for use 
c 

ql=qq(l,i,j,k) 

q2=qq(2,ij',k) 

q3=qq(3,i,j,k) 

q4=qq(4,i,j,k) 

q5=qq(5,i,j,k) 

c Compute pt/pO ratio 

ptp0(l)=gama*(gama-l)*q5+(gama-gama**2/2)/ql*(q2**2+q3**2+q4**2) 
c Compute T/Tt ratio and mach number 

tt0(l)=gama*(gama-l)/ql*(q5-04*(q2**2+q3**2+q4**2)/ql) 
m(l)=((l/u0(l)-1 )*2.0/(gama-1 ))**0.5 
c 

c Compute Theta and Non-D velocities theta, r and z 
c 

theta(l)=asin(y(ij,k)/(z(i,j,k)**2+y(i,j,k)**2)**0.5) 

vtnd(l)=-q4*sin(theta(l))+q3*cos(theta(l)) 

vmd(l)=q4*cos(theta(l))+q3*sin(theta(l)) 

vznd(l)=q2 

m(l)=(sqrt(vtnd(l)**2+vznd(l)**2))/ql 

degfl)=theta0)*57.296 


c Compute flow angle (degrees) using vtnd and vznd 

ang(l)=atan2(vtnd(l),vznd(l))*57.296 
15 continue 

write(7,300)um,deg( 1 ),ptpO( 1 ),m( 1 ),ang(l) 
write(8,300)um,deg(2),ptp0(2),m(2),ang(2) 
write(9,300)um,deg(3),ptp0(3),m(3),ang(3) 
write( 10,300)um ,deg(4),ptp0(4),m(4),ang(4) 
write(ll,300)um,deg(5),ptp0(5),m(5),ang(5) 

10 continue 


c Repeat above on other side of wake cut at imax 
c j=l not used (duplicate location across wake) 
c jmax not repeated (same location as j=31 at i=l) 
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jmaxm=jmax-l 
do 30 j=2,jmaxm,l 


do 25 1=1,5 
k=kk(l) 
ql=qq(l,i,j,k) 
q2=qq(2,i,jjc) 
q3=qq(3,i,j,k) 
q4=qq(4.i,j,k) 
q5=qq(5,ij,k) 

ptp0(l)=gama*(gama-l)*q5+(gama-gama**2/2)/ql*(q2**2+q3**2+q4**2) 
tl0(l)=gama*(gama-l)/ql*(q5-0.5*(q2**2+q3**2+q4**2)/ql) 
m(l)=(( l/ti0(l)- l)*2/(gama-1 ))* *0.5 
theia(l)=asin(y(i,j,k)/(z(i,j,k)**2+y(i,j,k)**2)**0.5) 
vtnd(l)=-q4*sin(theta(l))+q3*cos(iheia(l)) 
vmd(l)=q4*cos(theta(l))+q3*sin(theta(l)) 
vznd(l)=q2 

m(l)=(sqn(vtnd(l)**2+vznd(l)**2))/ql 
deg(l)=theta(l) *57.296 
ang(l)=atan2(vtnd(l),vznd(l))*57.296 
25 continue 

write(7,300)um,deg(l),ptp0(l),m(l),ang(l) 
write(8,300)um,deg(2),ptp0(2),m(2),ang(2) 
write(9,300)um,deg(3),ptp0(3),m(3),ang(3) 
write(10,300)um,deg(4),ptp0(4),m(4),ang(4) 
write(l 1,300)um,deg(5),ptp0(5),m(5),ang(5) 

30 continue 

: residual history output to unit 4 


write(4,310) l,(resd(l,l),l=l,5) 
do 40 j=2,nres 
it=10*0'-l) 

40 write(4,310)it,(resd(j,l),l=l,5) 

300 format(i5,5f8.3) 

310 format(i5,5( 1 x,e 10.3)) 

610 formate ***** warning *****',/, 

1 ' im, jm, km, read from input',3i5,’ do not 

2 ' im, jm, km, read from restart file',3i5) 


stop 

end 


match’,/. 
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APPENDIX H. MACH NUMBER AND FLOW ANGLE COMPARISONS 


Mach Number 

Mach Number Comparison 

1.78 mm Depth (89.8% Span) 



-8 - 

7 -6 -5 -4 -3 -2 -1 0 1 2 3 4 5 6 7 8 9 
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Figure HI. Mach Number Comparison At 1.78 mm Depth 



Position (Degrees) 


Figure H2. Mach Number Comparison At 0.42 mm Depth 
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Wake Flow Angle Comparison 

0.42 mm Depth (97.6% Span) 
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Figure H5. Flow Angle Comparison At 0.42 mm Depth 



Figure H6. Row Angle Comparison At 0.06 mm Depth 
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APPENDIX I. CONVERGENCE HISTORY 


RVC3D Residuals [0.5070 Pressure Ratio] 

1.00E-04 
■2 1.00E-05 
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J 1 1.00E-07 
1.00E-08 

1 500 1000 1500 2000 2500 3000 

Number of Iterations 

Figure II. 0.5070 Pressure Ratio Convergence History 



RVC3D Residuals [0.6041 Pressure Ratio] 



Number oflterations 


Figure 12. 0.6041 Pressure Ratio Convergence History 
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RVC3D Residuals [0.9054 Pressure Ratio] 



Number of Iterations 


Figure 15. 0.9054 Pressure Ratio Convergence History 
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APPENDIX J. MULTIPLE GRID AND SOLUTION PLOTTING 


Program mgrid.f read an x, y, z formatted grid file (fort.l) and converted it into a 
multiple grid file format (fort.21) for use with PLOT3D. 


Program mgrid.f 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c Program to read a single PLOT3D file and convert 
c it to a mgrid file format 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
real x(250,65,65).y(250,65,65),z(250.65,65) 
real yy(250,65,65),zz(250,65,65) 
read( 1 ,*)idim,jdim,kdim 

read(l,*)(((x(i,j,k),i=l,idim)j=l,jdim),k=l,kdim), 

# (((y(i,j,k),i=l,idim),j=l,jdim),k=l,kdim), 

# (((z(ij,k),i=l,idim).j=l,jdim),k=l,kdim) 
write(*,*)'Done reading fort.l' 

c calculate pitch 
jm=jdim 

"r2=y(ljm.l) + *2+z(l,jm,l)**2 
sinp=(y(l ,jm,l )*z(im,jm,l )-y(im jm, l)*z(l jm.l ))/r2 
pitch=asin(sinp) 
c 

do i=l,idim 
doj=l jdim 
do k=l,kdim 

r=(z(i,j,k)**2+y(i,j,k)**2)**0.5 

lheta=asin(y(i,j.k)/r) 

theta=theta+pitch 

yyCi.j ,k)=r*sin(theta) 

zz(i ,j ,k)=r*cos(theta) 

enddo 

enddo 

write(21)2 

write(21 )(idim,jdim,kd im,igrid= 1,2) 
c do 10 igrid=l,2 

write(21)(((x(ij,k)4=l,idim)j=l jdim),k=l,kdim), 

# (((y(i,j,k),i=l,idim).j=l ,jdim),k=l,kdim), 

# (((z(i j,k),i=l ,idim)u=l jdim),k= 1 ,kdim) 
c 10 continue 

write(21)(((x(ij.k),i=l,idim),j=l,jdim),k=l,kdim), 

# (((yy(i,j,k),i=l ,idim),j=l ,jdim).k= 1 ,kdim), 

# (((zz(i,j,k),i= 1 .idirn) ,j=1 jdim),k= 1 ,kdim) 
stop 

end 
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Program mq.f read a solution file (fort.3) and converted it into a multiple solution file 
format (fort.23) for use with PLOT3D. 


Program mq.f 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

c Program to read a single PLOT3D file and convert 
c it to a Q_mgrid file format 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
real qq(250,65,65.5) 
read(3,*)idim,jdim,kdim 
read(3 ,*)fsmach,alpha,re,time 

read(3,*)((((qq(i,j,k,nx),i= 1 ,idim) j=l jdim),k=l ,kdim),nx= 1.3) 
write(23)2 

write(23)(idim.jdim,kdim,igrid=l,2) 
do 10 igrid=l,2 
write(23)fsmach,alpha,re,time 

write(23)((((qq(i,j,kjtx),i=l,idim)j=l,jdim),k=l,kdim),nx=l,5) 

10 continue 
stop 
end 
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